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The  Molecular  Modelling  of  Silane  Primers  for  Adhesive 
Joints  and  Composite  Materials 


1.  Introduction 

Recently,  experimental  work  has  been  performed  in  the  Department  of  Mechanical 
Engineering  at  Imperial  College  of  Science,  Technology  and  Medicine  by  N.G.  Cave 
and  AJ.  Kinloch^*^  investigating  the  use  of  self-assembling  monolayer  silane  films  as 
adhesion  promoters.  This  was  followed  up  by  computational  work  of  M.S.  Sennett'^^ 

et  al.  and  S.A.  Ojo‘^^  examining  the  interaction  of  hydrolysed  y-glcidoxypropyltrimeth- 

oxysilane  (GPMS)  and  nine  variants  thereof,  on  four  metal  oxide  surfaces.  We  intend 
to  continue  on  from  this  work,  and  initially  attempt  to  use  computer  modelling  methods 
developed  in  references  two  and  three  to  reproduce  the  experimental  results  in  reference 
one. 

2.  Aims 

Hie  main  aim  of  this  wolk  is  to  improve  the  understanding  of  the  molecular-level 
mechanism  of  adhesion  of  organosilanes  on  metal  surfaces  for  the  design  of  improved 
adhesion  promoters/primers.  Our  approach  will  be  to  carry  out  molecular  simulations 
to  determine  binding  energies  of  the  various  molecular  structures,  and  the  influence  of 
surface  contaminants,  on  the  interaction  of  organosilanes  with  metal  oxides.  We  will  be 
using  a  variety  of  computer  simulation  methods  incorporated  in  the  Cerius^  software 
package^'**,  which  are  discussed  below.  We  will  correlate  the  models  with  experimental 
data  in  order  to  validate  the  model. 


3.  Computer  Simulation  Methods 
3.1  Minimisation 

Minimisation  involves  a  mathematical  algorithm  to  find  the  lowest  energy  of  a  given 
structure.  During  a  minimisation,  the  atomic  co-ordinates  are  adjusted  in  order  to 
reduce  the  molecular  energy.  The  energy  of  the  system  is  calculated  from  a  series  of 
force  field  terms  that  describe  the  interactions  between  the  atoms.  The  interactions 
include  bond  terms,  angle  terms,  torsion  terms,  non-bonded  interactions,  and  charge 
terms.  In  our  calculations,  the  metal  oxide  surface  will  remain  fixed.  Initially  we 
intend  to  take  our  silane  molecules,  and  minimise  them  onto  the  surface,  then  examine 
the  interactions. 


3.2  Molecular  Dynamics  (MD)  Calculations 

Molecular  Dynamics  is  an  adaptation  of  the  minimisation  process.  A  minimisation  will 
take  the  molecule  to  the  nearest  lowest  energy,  which  is  not  necessarily  the  global 
minimum.  In  order  to  overcome  this,  a  molecular  dynamics  calculation  can  be 
performed.  This  places  an  initial  velocity  on  each  of  the  atoms,  which  is  related  to  a 
temperature  factor.  The  atoms  are  allowed  to  move  for  a  given  time  step,  within  the 
limits  allowed  by  the  force  field  terms,  then  the  energy  of  the  new  structure  is  calculated 
using  the  force  field  terms,  as  in  a  minimisation  calculation.  By  allowing  the  structure 
to  be  ‘heated’,  it  can  overcome  the  energy  barriers  between  localised  minimum,  and 
allow  a  number  of  low  energy  positions  to  be  observed.  The  structure  after  each  time 
step  is  saved,  allowing  a  trajectory  file  to  be  saved,  which  can  be  analysed  at  a  later 
time.  It  is  intended  to  run  molecular  dynamics  calculations  on  the  various  molecules,  to 
determine  minimum  energy  positions  and  likely  conformations  and  interactions. 

3.3  Sorption  and  Minimisation 

The  sorption  process  allows  molecules  to  be  placed  at  random  orientations  above  the 
metal  oxide  surface.  An  algorithm  selects  one  of  four  methods  to  generate  new 
positions 

•  Creation  of  a  new  position,  according  to  a  collision  check,  then  an  energy  selection 
criteria  based  on  probability. 

•  Destruction  of  a  molecule  -  a  molecide  is  removed  according  to  the  same 
energy/probability  criteria  as  before. 

•  Rotation  of  an  existing  molecule  by  a  random  amount  within  a  pre-defined  limit. 

•  Translation  of  an  existing  molecule,  again  by  a  random  amount  within  a  pre-defined 
limit. 

The  energy  is  calculated  for  each  molecule,  and  the  position  recorded,  as  in  the  previous 
dynamics  calculations.  This  method  allows  us  to  build  up  a  picture  of  many  molecules 
on  a  surface,  and  how  they  are  likely  to  pack  together.  It  can  also  be  used  to  generate 
up  starting  position  for  dynamics  runs.  Analysis  of  the  results  will  allow  us  to  build  up 
a  map  of  how  the  interaction  energy  changes  with  distance  from  the  surface. 

4.  Conclusions 

Initial  work  has  concentrated  on  modelling  basic  straight  chain  alcohols  from 
propan-l-ol  to  dodecan-l-ol.  This  is  to  enable  the  validation  of  a  force  field  developed 
by  Andreas  Bick  at  Biosym/MSf*’,  with  respect  to  the  hydrogen  bonding  that  is 
expected.  Once  this  has  been  validated,  we  will  proceed  on  to  examine  the  silanes. 
Once  we  have  validated  the  force  field  against  experimental  values,  we  will  then 


increase  the  complexity  of  onr  calculations  to  include  water  molecules.  Once 
successful,  we  wiU  then  go  on  to  develop  suggestions  of  new.  improved  pnmer  types, 

based  on  our  findings. 
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